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AN INTEGRAL REPRESENTATION OF THE GENERALIZED 
EULER-MASCHERONI CONSTANTS 


The generalized Euler-Mascheroni constants are defined by 
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and are coefficients of the Laurent expansion of the Riemann Zeta function f(z) about z = 1 
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These were first defined by Stieltjes in 1885, discussed by Stieltjes and Hermite [1], and have been 
periodically reinvented over the years [2,3,4]. Computation of these numbers has attracted several 
writers, among them: Jensen [5], n=l to 9; Gram [6], n=l to 16; Liang and Todd [7], n=l to 20; 
and Ainsworth and Howell [8], n=l to 31. Reference 7 uses Euler-Maclaurin summation with an 
extended precision software package but gives no error bounds. Reference 8 extends the list by another 
1 1 and does almost as well in terms of accuracy as Reference 7 using standard precision and some 
mathematical analysis and does include error bounds. 

However it seems that 35 is about the limit that one can reach using Euler-Maclaurin summation 
of the series as given in equation (1) so we shall express this series as an integral which is easy to 
numerically evaluate. This integral will allow us to compute the first 2000 Euler-Mascheroni constants. 

Consider the integral 
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taken about the contour C in Figure la (Fig. lb depicts the magnitude of the integrand for n=4) where 
0 < y < R, 1 < x < N + 1/2. Further, we take 
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Now consider the first part of I2, I3, and I4 (the second, third, and fourth integrals). We have: 
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and the sum 


ln n+1 (N+l/2) 
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may be replaced by 

ln n+1 (N+1 ) 
n+1 

when we let N approach infinity. 

We now have, upon slight rearrangement, 
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Because of the e 2,rR factor in the denominator, the last two integrals obviously tend to zero as 
R approaches infinity. Also, the first two integrals are conjugates so that 
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which tends to zero as N approaches infinity. Therefore, 
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As a check, we make the substitution 
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in the series portion of equation (2) to obtain: 
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Thus, we obtain for the Riemann zeta function 
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in agreement with the formula derived by the Danish mathematician J. L. W. V. Jensen in 1893 [9] 


COMPUTATION 


Using the new representation 
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Now, for n > 150, the integrand becomes too large for computational purposes. We thus introduce a 
scale factor 0 n which can be chosen to keep the integrand well behaved. Thus, we obtain 



After the integral is evaluated, it can be easily multiplied by the scale factor to give the desired result. 

The Gauss numerical integration formula was used and if the value was unchanged by taking 
n=N and n=2N, we assumed the integral to be correct as far as the agreement extended. The integrand 
is much too complex to use the standard remainder terms and the bounds that we were able to find in 
Reference 8 were very much too conservative. A further check is that the results were in agreement 
with those obtained when rectangular integration was used. 

Table 1 gives 7 n , n = 1 (1) 200, and assorted values to n = 2000. 

We now take n in 7 n to be a non-negative continuous variable and obtain the graph (Fig. 2) 
plus a listing of the first 13 zeros (Table 2). Conceivably the integral of 7 n or some function of its 
roots might lead to some esoteric information regarding the Riemann zeta function. 
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TABLE 1 . GENERALIZED EULER-MASCHERONI CONSTANTS 


n 


n 

T'n 

n 

0 

5.7721 56649E-1 

40 

2.487215594E-1 

80 

1 

-7.281584548E-2 

41 

-7.195748469E-1 

81 

2 

-9.690363 193E-3 

42 

-2.638794927E0 

82 

3 

2.053834420E-3 

43 

-5.2649303 12E0 

83 

4 

2.325370065E-3 

44 

-7.188745890E0 

84 

5 

7.933238173E-4 

45 

-5.072344590E0 

85 

6 

-2.387693454E-4 

46 

6.6099 15609E0 

86 

7 

-5.272895671 E-4 

47 

3.403977498E1 

87 

8 

-3.521233538E-4 

48 

7.868247976E1 

88 

9 

-3.439477443E-5 

49 

1.258443876E2 

89 

10 

2.053328149E-4 

50 

1.268236027E2 

90 

1 1 

2.701 844395E-4 

51 

-1.919691 183E1 

91 

12 

1.67272912 1 E-4 

52 

-4.63 1889230E2 

92 

13 

-2.746380660E-5 

53 

-1.340659144E3 

93 

14 

-2.092092621 E-4 

54 

-2.57245474 1E3 

94 

15 

-2.834686553E-4 

55 

-3.457141209E3 

95 

16 

-1.996968583E-4 

56 

-2.055275816E3 

96 

17 

2.6277037 1 1 E-5 

57 

5.372282214E3 

97 

18 

3.073684081 E-4 

58 

2.401938938E4 

98 

19 

5.036054530E-4 

59 

5.742431929E4 

99 

20 

4.66343 56 15 E-4 

60 

9.854325459E4 

100 

21 

1.044377698E-4 

61 

1.1 16709578E5 

101 

22 

-5.41 5995822E-4 

62 

5.333665226E3 

102 

23 

- 1.243 96 2090 E-3 

63 

-3.909726873E5 

103 

24 

-1.58851 1279E-3 

64 

-1.3031 807 13E6 

104 

25 

-1.074591953E-3 

65 

-2.845076553E6 

105 

26 

6.568035 186E-4 

66 

-4.540526610E6 

106 

27 

3.477836914E-3 

67 

-4.341 905 139E6 

107 

28 

6.400068532E-3 

68 

2.871566946E6 

108 

29 

7.3711 5 1770E-3 

69 

2.660490855E7 

109 

30 

3.557728856E-3 

70 

7.9321 663 12E7 

110 

31 

-7.5 13325998E-3 

71 

1.6621 5 1340E8 

111 

32 

-2.57037291 1 E-2 

72 

2.55 1532583E8 

112 

33 

-4.51 06734 11E-2 

73 

2.126556317E8 

113 

34 

-5.1 12692802E-2 

74 

-2.987670894E8 

114 

35 

-2.037304360E-2 

75 

-1.919487428E9 

115 

36 

7.2482 15882E-2 

76 

-5.515574258E9 

116 

37 

2.360263823E- 1 

77 

-1.148345099E10 

117 

38 

4.2896344641’’- 1 

78 

-1 .75701 5228E10 

118 

39 

5.1 792 1 84271-:-! 

79 

- 1 .396 1021 46E 1 0 

119 


2.5 16344101E10 
1.510585108E1 1 
4.379044312E1 1 
9.317068469E1 1 
1 .47209981 9E 1 2 

1.259044968E12 

-1.958810225E12 

-1.295154550E13 

-3.929715388E13 

-8.753044740E13 

-1.471610494E14 

-1.477977001E14 

1.294632141E14 

1.188562976E15 

3.920684627E15 

9.341670850E15 
1.707524745E16 
2.071798354E16 
-2.854307850E1 5 
-1.125 848 1 08E 1 7 

-4.253401 572E17 
-1.107785955E18 
-2.232843819E18 
-3.249016082E18 
-1.590464731E18 

1 .01 782 1085E19 
4.868854920E19 
1.424392927E20 
3.200736262E20 
5.487745762E20 

5.449976660E20 
-6.7290 19783E20 
-5.61 1293844E21 
-1.930664135E22 
-4. 902622623 E22 

-9.759599993E22 
-1.373492481 E23 
-3.817890627E22 
5.8968478471-23 
2.656742888E24 


8 



TABLE 1. (Concluded) 


n ? n 

120 7.81 8403277E24 

121 1 .79646 1571 E25 

122 3.181276594E25 

123 3.336371622E25 

124 -3.782268183E25 

125 -3.4755161 17E26 

126 -1.257375875E27 

127 -3.362216797E27 

128 -7.131871 952E27 

129 -1.1 1532378 1E28 

130 -6.637629576E27 

131 3.605865029E28 

132 1.936502639E29 

133 6.244838068E29 

134 1.564173210E30 

135 3.101716787E30 

136 4.22670 1286E30 

137 -7.995222356E27 

138 -2.523130101 E3 1 

139 -1.103550302E32 

140 -3.323756372E32 

141 -7.95305091 1E32 

142 -1.490838047E33 

143 -1.747042301E33 

144 1.370337533E33 

145 1.679138192E34 

146 6.623981071E34 

147 1.914935612E35 

148 4.445803239E35 

149 7.981621278E35 

150 8.028853731E35 

151 -1.396229273E36 

152 -1.134586466E37 

153 -4.241080464E37 

154 -1.199627658E38 


n Tn 

155 -2.740591577E38 

156 -4.790923 159E38 

157 -4.24774747 1E38 

158 1.143679930E39 

159 7.9707423 18E39 

160 2.910661837E40 

161 8.179814602E40 

162 1 .86381 9327E4 1 

163 3.236126153 E4 1 

164 2.703985771E41 

165 -8.81 1044847E41 

166 -5.871467633E42 

167 -2. 142378 188E43 

168 -6.061386654E43 

169 -1.395436964E44 

170 -2.457784157E44 

171 -2.124640956E44 

172 6.605630678E44 

173 4.538504341E45 

174 1.687141829E46 

175 4.863746403E46 

176 1. 144982704E47 

177 2.086041054E47 

178 2.021973061E47 

179 -4.769958441E47 

180 -3.664699738E48 

181 -1.414931352E49 

182 -4.203887079E49 

183 -1.022891247E50 

184 -1.958109259E50 

185 -2.218374141E50 

186 3.1 12089130E50 

187 3.061460410E51 

188 1 .25525342 1E52 

189 3.887946022E52 




190 

9.870435858E52 

191 

2.0071 16958E53 

192 

2.676425 145E53 

193 

-1.391 163860E53 

194 

-2.601448677E54 

195 

-1.167447626E55 

196 

-3.81 7302910E55 

197 

-1.019713344E56 

198 

-2.217662282E56 

199 

-3.4420708 18E56 

200 

-6.974649720E55 

250 

3.059212855E79 

300 

-5.556728220E102 

350 

-1.785446107E127 

400 

-1.761642187E152 

500 

-1.165505188E204 

750 

-2.1645871E341 

1000 

-1.5709538E486 

1500 

-3.1 1 70 1 39E79 1 

2000 

2.68024593E1 109 
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TABLE 2. ROOTS p [ LESS THAN 60 OF 7 , WHERE I < 10 


2.623384858 

5.684688574 

9.114834513 

12.86713391 

16.90281884 

21.18896 

25.6981707 

30.40797324 

35.29985207 

40.35837636 

45.57051042 

50.9251026 

56.41250687 



Figure 2. The Euler-Maseheroni function. 
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